home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
IRIX 6.2 Development Libraries
/
SGI IRIX 6.2 Development Libraries.iso
/
dist
/
complib.idb
/
usr
/
share
/
catman
/
p_man
/
cat3
/
complib
/
DPBSL.z
/
DPBSL
Wrap
Text File
|
1996-03-14
|
2KB
|
67 lines
DDDDPPPPBBBBSSSSLLLL((((3333FFFF)))) DDDDPPPPBBBBSSSSLLLL((((3333FFFF))))
NNNNAAAAMMMMEEEE
DPBSL - DPBSL solves the double precision symmetric positive definite
band system A*X = B using the factors computed by DPBCO or DPBFA.
SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
SUBROUTINE DPBSL(ABD,LDA,N,M,B)
DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
On Entry
AAAABBBBDDDD DOUBLE PRECISION(LDA, N)
the output from DPBCO or DPBFA.
LLLLDDDDAAAA INTEGER
the leading dimension of the array ABD .
NNNN INTEGER
the order of the matrix A .
MMMM INTEGER
the number of diagonals above the main diagonal.
BBBB DOUBLE PRECISION(N)
the right hand side vector. On Return
BBBB the solution vector X . Error Condition
AAAA division by zero will occur if the input factor contains a zero on the
diagonal. Technically this indicates singularity, but it is usually
caused by improper subroutine arguments. It will not occur if the
subroutines are called correctly, and INFO .EQ. 0 . To compute
INVERSE(A) * C where C is a matrix with P columns
CALL DPBCO(ABD,LDA,N,RCOND,Z,INFO)
IF (RCOND is too small .OR. INFO .NE. 0) GO TO ...
DO 10 J = 1, P
CALL DPBSL(ABD,LDA,N,C(1,J)) 10 CONTINUE LINPACK. This version dated
08/14/78 . Cleve Moler, University of New Mexico, Argonne National Lab.
Subroutines and Functions BLAS DAXPY,DDOT Fortran MIN0
PPPPaaaaggggeeee 1111